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In a subsequent analysis a coding error was discovered in the treatment of the short range correlations. Correcting 
the error resulted in an increase of the neutrinoless double beta decay matrix elements. Here we provide the most 
relevant correct numerical results. The corrected Table (replacing Table 1 of the paper) and the corrected figure 
(replacing Fig. 2 of the paper) are shown. 

While the matrix elements are now larger, our basic claim that the chosen way of adjusting the interaction strength 
makes the matrix elements essentially independent on the size of the single particle basis, on the parametrization of 
the G-matrix, whether QRPA or RQRPA is used (although, as seen, QRPA results in ~ 10% larger matrix elements 
than RQRPA), and whether gA is quenched or not remains true. 

TABLE I: Averaged 0v/3f3 nuclear matrix elements (M l0v ) and their variance a (in parentheses) evaluated in the RQRPA and 
QRPA. In column 6 the variance e exp . of the O^/3/3-decay matrix element due to uncertainties in the measured 2^/3/3-decay 
half-live T^~ elcv is given. and gA denote the 2i//3/3-decay nuclear matrix element deduced from Ty 2 exp and axial- vector 

coupling constant, respectively. In column 7 the Ov (3/3 half- lives evaluated with the RQRPA average nuclear matrix element 
and for assumed {mpp) — 50 meV are shown. For 136 Xe there are four entries; the upper two use the upper limit of the 2v 
matrix element while the lower two use the ultimate limit, vanishing 2v matrix element. 150 Nd is included for illustration. It 
is treated as a spherical nucleus; deformation will undoubtedly modify its fV matrix element. 
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FIG. 1: Average nuclear matrix elements (M 1 v ) and their variance (including the uncertainty coming from the experimental 
error in M 2v ) for both methods and for all considered nuclei. For 136 Xe the error bars encompass the whole interval related to 
the unknown rate of the 2v/3/3 decay. 



